Small— Angle Bhabha Scattering at LEPl. 
Analytical Results for Wide— Narrow Angular 

Acceptance. 

A.B. Arbuzov% G.I. Gach^, V.Yu. Gontchar^, E.A. Kuraev^, 
N.P. Merenkov^ and L. Trentadue*^ 



' ^Bogoliuhov Laboratory of Theoretical Physics, JINR, Duhna, 141980, Russia 

0^ ' ^Institute of Physics and Technology, Kharkov, 310108, Ukraine 

^ \ ^ Dipartimento di Fisica, Universitd di Parma and INFN, Gruppo Collegato di 
Q ' Parma, 43100 Parma, Italy 

^ : 



■ Abstract 
> , 

i Analytical formulae for small-angle Bhabha scattering cross-section at LEPl are 

CN ' given for the case of wide-narrow angular acceptance. Inclusive and calorimeter 

■ event selections are considered. Numerical results are presented. 

Q^ ■ PACS 12.20.-m Quantum electrodynamics, 12.20.Ds Specific calculations 
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^ i Introduction 



The small-angle Bhabha scattering (SABS) process is used to measure luminosity of electron- 
positron colliders. At LEPl an experimental accuracy on the luminosity better then one per 
mille has been reached [1]. To estimate the total accuracy a systematic theoretical error must 
be added. That is why in recent years a considerable attention has been devoted to theoretical 
investigation of SABS cross-section [2-11]. 

The theoretical calculation of SABS cross-section at LEPl has to cope with two problems. 
The first one is the description of experimental restrictions in terms of final particles phase 
space used for event selection. The second one concerns the computation of matrix element 
squared with the required accuracy. There are two approaches for the theoretical study of 
SABS at LEPl: the one basing on Monte Carlo (MC) programs [2-6] and the other using 
semi-analytical calculations [7-11]. 
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The advantage of MC method is the possibihty to model different types of detectors and 
event selections [2]. But at this approach some problems with exact matrix element squared 
exist. Contrary, the advantage of the analytical approach is the possibility to use exact 
matrix element squared and its defect is a low mobility relative the change of an experimental 
conditions for event selection. Nevertheless, the analytical calculations are of great importance 
because they allow to check numerous MC calculations for different types of ideal detectors. 

In this letter we list some analytical results for SABS cross-section at LEPl suitable for 
inclusive (labeled in [2] as BAREl) and calorimeter (CALOl and CAL02) event selections 
in the case of asymmetrical wide-narrow circular detectors. We give analytical formulae for 
the full first order correction to the cross-section as well as for leading second and third order 
ones. Our numerical estimations include also next-to-leading second order contribution in 
the case BAREl. 



1 First order corrections 



Let us consider at first BAREl event selection. We introduce the dimensionless quantity 
S = Q^o'obs/(47ra^), where Qf = e^Of {e is the beam energy and 6i is the minimal scattering 
angle for wide circular detector) and Uobs is an experimentally observable (by means of ideal 
detectors) cross-section. Then the first order QED correction can be written as follows: 
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where 



K{x, z- p4, P2) = + Ag) + A42L1 + AgL2 + (ef - e(^))L3 



+ (64-62)^4, 

K{x, P3, 1) = + + ^1 + + ef L3 , (2) 

1 + 

and L — ln{zQl/ml) is the large logarithm; parameter Xc puts the restriction on the energies 
of the final electron and positron: 8182 > XcS^. In Eqs.(l) and (2) we used the following 
notation for 0-functions and Lj: 
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The quantities Lj can be obtained from Lj by the substitution Pi ^ Pa, p2 ^ I- Here we use 
the same notations as in [10]. Note only that for wide-narrow angular acceptance 



P3 > P4 > P2 > Pi = 1 , Pi^T^ (3) 

where 9i are the hmiting angles of the circular detectors (see Section 3). Function P'^^\x) 
defines the iterative form of the non-singlct electron structure function (sec for example [6]). 
The first (second) line in the r.h.s. of Eq.(l) is the contribution due to photon emission by 
positron (electron) . The terms accompanied with x-dependent 0-functions under the integral 
sign correspond to the initial-state correction while the rest belongs to the final-state one. 

The CALOl cluster is the cone with angular radius 6 = 0.01 around the direction of the final 
electron (or positron) momentum. If photon belongs to the cluster the whole cluster energy is 
measured, and electron may have any possible energy. Therefore, the limits of x-integration 
for Eobs extend from to 1 here. If photon escapes the cluster the event looks the same as 
for BAREl. The above restrictions on x-integration limits can be written symbolically as 
follows: 

1 Xc 1 Xc ^ 

Jdx + e{5-\r\)Jdx ^ J dx - e{\r\-S) J dx, r = --^, (4) 

Xc LO El 



where k (q^) and co (si) are the transverse momentum and the energy of the hard photon 
(electron). 

As we saw on the example of BAREl event selection it is necessary to distinguish the con- 
tributions into El due to electron and positron radiation: 

El = + . (5) 



According to (4) we have 

Et = E, + E/ + E; + E} , E^ = Ei + E^ + E; + E; , (6) 

where index i{f) labels the initial (final) state and index c points on a cluster-form depen- 
dence. 
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Quantities Sj and coincide with the corresponding initial-state correction for BAREl (see 
comments on Eq.(l)). For T,f and we can use the form of differential cross-section suitable 
for inclusive event selection with extended x-integration limits: 
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In order to find the additional contributions into Si which depend on the cluster form it is 
enough to use the simplified differential cross-section of single photon radiation, neglecting 
electron mass, and taking into account the restrictions |r| < 6 (for the initial state) and 
|r| > 5 (for the final state). The contribution due to photon emission by the initial-state 
electron can be written as follows: 





where defines the integration limits over z (in the square brackets) and over zi (in the 
parenthesis) : 

[a^, al]{x'^z+,x'^) + [6^, a'^]{x'^z+,x^Z-) + [bl, b^]{x^pl,x^Z-), ao = P2, &o = P4, 
a = max(p2, 1 + A(l — x)), 6 = min(p4, p3 — A(l — x)), z± — {y/z ± X{1 — x)'^) , 

and function $ is defined below: 
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The cluster-dependent contribution due to photon emission by the final-state electron reads 
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To obtain it is enough to substitute in the expression for ^' the parameters a, b, ao and 
60 by d, b, So and 60 respectively, 



d = p2 + A(l — x), ao = max(l, p2 — A(l — x)), 
b = p4- X{1- x), bo = min(p4 + A(l - x),p3). 

Finally, the cluster-dependent contribution due to photon emission by the final-state positron 
can be written as follows: 

1 

+ jff{a,b,ao,bo ^ d,b,do,bo; Ps, Pi ^ P4, P2) ■ (10) 



The CAL02 event selection differs from the CALOl one by the form of the cluster (see [2]). 
Only cluster-dependent contributions into Si will be changed here. Analytical formulae are 
very cumbersome, and we give the result for symmetrical wide-wide case only (S''' = S^) but 
our numerical calculations include wide-narrow angular acceptance too. 
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We present also the corresponding formula for photon emission by the final-state electron: 
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Fi = arctan ^^py| , -F2 = arctan(C), F3 = arctan((5^^^), C = T/cot ' 



Angle $ and quantity A, which enter into Eq.(12), define the form and the size of CAL02 
cluster. Namely [2] 



Stt - ^0 . _ 0.051 



Finally, functions J(±) and z^"^^ are defined as follows: 



Ji 



zp - x\l - xfX^ sin2$ ± (1 - x)A(l - 2xsin2 



,2 ^ (±) 



13 = 1- Ax{l - x) sin" -, zl^> = {pi ± (1 - x)\f - Ax{l - x)pi{pi ± A) sin" - 



2 Second and third order corrections 



In this Section we give the analytical form of the leading second and third order corrections 
to SABS cross-section suitable for both, inclusive and calorimeter, event selections. The 
contribution connected with pair production in the singlet channel is negligible for LEPl 
conditions, while the one in the non-singlet channel can be taken into account by means of 
effective QED couphng [6] . Therefore, we will consider here the photonic corrections only. 

The second order correction can be presented in the form 



(13) 
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The first term in r.h.s. of Eq.(13) is responsible for two-photon (real and virtual) emission by 
the electron. The second one describes two-photon emission by the positron. And the third 
one considers the situation when both the electron and the positron radiate. 

The leading contributions in the case of inclusive event selection read 
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where 



1 1 

p(2)(^) = I ^p(l)(t)p(l) (^1^ , I p(^)(x)dx = 0. 

X 

The contribution of initial-(final-)state radiation for J^'^'y^ and S:^^ are accompanied with 
x-dependent (^-independent) ©-functions and xi, a;2-dependent (xi, X2-independent) ones 
for E^^. The terms with additional integration over t- variable describe the simultaneous 
radiation of one photon from the initial state and the other from final state (initial-final- 
state radiation). 

In the case of calorimeter event selection we have to take in the r.h.s. of Eqs.(14) the terms 
corresponding to initial-state radiation only. The elimination of final-state one exhibits itself 
by means the last equality. As concerns the contribution due to initial-final-state radiation 
it may be understood as follows. 

In fact i-variable in Eqs.(14) is the energy fraction carried out by both the final-state radiated 
photon and the final-state electron (or positron). Just this value defines the cluster energy for 
calorimeter event selection. The x-variablc which is the energy fraction of the final electron 
by definition can change here from up to t. That is why initial-final-state radiation of the 
electron for calorimeter event selection will be proportional to 

/dt/^p«(t)p«(^)e?) = /dtp«(t)0j)/p«(i/)d|/ =0. (15) 

Xc Xc 



The same is valid of course for the corresponding part of positron radiation. 
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Thus, we see that it is enough to have only single final-state radiated photon to eliminate the 
leading contribution due to initial-final-state radiation. This conclusion refiects the essence 
of a reduced Lee-Nauenberg theorem [12] and is valid for all higher order corrections. 

The leading third order correction reads 



2 oo 1 1 

Xc Xc/x 



1 rdt 

+ 



X 

Z, = i[(A3i + Aif)(A42 + a£^)) + (A31 + Air))(A42 + Ag)]p(^)(a;)p(^)(xi) 
+ pW(x) / y [Ag(A42 + Ag) + Ag(A3i + Ag)]P«(t)pW (^) , 

XI 

p(3)=|^p(i)(t)p(2) (^^y 



For calorimeter event selection it is needed to take 

Z, = ^(AgAai + Ag) A42)P(^)(^), = ^(A^? A&) + A^^,^^A^S)P^'\x)P^'\x,). 

Note that in this case the leading second and third order contributions have a universal 
character and do not depend on cluster form. Thus, they are suitable for both, CALOl and 
CAL02. 

Quantity Zi describes the situation when only one fermion (electron or positron) radiate 
(one-side emission), while Z2 is responsible for simultaneous radiation of the electron and 
the positron (opposite-side emission). 

The formulae for leading second and third order contributions are written in the form with 
different ©-functions under integral sign. One can eliminate these ©-functions using such 
kind relations as, for example 

J eSt^e^s'^dz dx dxi ^ J dz J dx j dxi. (17) 

^cP3 XcPzl\/z Xc/x 

It is needed to keep in mind only that every integral has to be equal to zero if the lower limit 
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of ^-integration becomes more than the upper one. The last statement is vahd for the first 
order correction too. 



3 Numerical results 

In our calculations we restrict ourselves with pure QED corrections supposing Z-exchange, 
vacuum polarization and up-down interference are switched off. As shown in papers by 
W. Beenakker and B. Pietrzyk [7], a sufficiently accurate luminosity determination requires 
the full Born plus complete order a corrected cross-section. Nevertheless our numerical re- 
sults can be used for comparisons and cross-checks with numerous Monte Carlo and semi- 
analytical computations [2]. 

We performed the numerical calculations for the beam energy — 46.15 GeV and the 

following sets of limiting angles of circular detectors: 



(limits of integration are (P4,pi) for NN and WN angular acceptances and {p\, 1) for WW 
one) equals to 

175.922 nb BAREl, CALOl WW, 

139.971 nb BAREl, CALOl NN; CAL02 WW , 

103.299 nb CAL02 NN . 

Formula (18) takes into account the contributions of scattering diagram and interference 
of scattering and annihilation ones. The contribution of pure annihilation diagram is pro- 
portional to 6\. It is neghgible even at the born level. When calculating the QED correc- 
tions to cross-section (18) we ignore systematically the terms proportional 9\. Terms of 
this kind have double logarithmic asymptotic behavior [13] and are equal parametrically to 
{oi\t\) W(|i|/s)/(7rs), what is about 0.1 per mille as compared with unit for LEPl conditions. 

The results of our calculations of QED correction with the switched off vacuum polarization 
are shown in the Table 1. The centre-of-mass energy is ^/s = 92.3 GeV. The second order 



i) BAREl, CALOl: = 0.024, 9^ = 0.058, 

ii) CAL02: ^1 = 0.024 + /i, 9^^ 0.058 - h, 



9i + h, ^4 = ^3-/1, 

9i + h, ^4 = ^3 - 3/i . 



0.017 
8 



The Born cross-section 




(18) 
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Table 1 



The results of our analytical calculations for the SABS cross-section 







BAREl 






CALOl 






CAL02 






WW 


NN 


WN 


WW 


NN 


WN 


WW 


NN 


WN 


a (nb) with 0{a^) corrections 


0.1 


166.008 


130.813 


134.504 


166.285 


131.032 


134.270 


130.997 


94.666 


98.354 


0.3 


164.702 


129.797 


133.416 


166.006 


130.833 


134.036 


130.705 


94.491 


98.127 


0.5 


162.203 


128.001 


131.428 


165.244 


130.416 


133.466 


130.141 


94.177 


97.720 


0.7 


155.390 


122.922 


125.809 


161.749 


128.044 


130.542 


127.491 


92.981 


95.874 


0.9 


134.334 


106.478 


107.945 


149.866 


118.822 


120.038 


117.491 


86.303 


87.696 






a 


(nb) with 0{a^) and 0{a^) photonic corrections 






0.1 


166.958 


131.674 


134.808 


167.073 


131.740 


134.572 


131.705 


95.334 


98.609 


0.3 


165.447 


130.534 


133.583 


166.686 


131.467 


134.231 


131.339 


95.118 


98.314 


0.5 


162.574 


128.474 


131.127 


165.718 


130.903 


133.471 


130.628 


94.731 


97.793 


0.7 


155.597 


123.206 


125.255 


162.042 


128.361 


130.378 


127.808 


93.377 


95.782 


0.9 


137.153 


108.820 


109.677 


150.732 


119.560 


120.411 


118.229 


86.931 


87.961 


absolute values of the 0(0^) pair production correction (nb) 


0.1 


0.007 


-0.004 


0.015 


-0.046 


-0.045 


-0.024 


-0.045 


-0.047 


-0.024 


0.3 


-0.033 


-0.033 


-0.020 


-0.046 


-0.045 


-0.024 


-0.045 


-0.047 


-0.024 


0.5 


-0.058 


-0.050 


-0.041 


-0.048 


-0.046 


-0.025 


-0.046 


-0.047 


-0.024 


0.7 


-0.090 


-0.074 


-0.069 


-0.069 


-0.059 


-0.042 


-0.059 


-0.051 


-0.036 


0.9 


-0.142 


-0.115 


-0.115 


-0.137 


-0.111 


-0.102 


-0.111 


-0.085 


-0.075 



absolute value of the O(a^) leading correction (nb) 



0.1 


-0.055 


-0.047 


-0.006 


-0.041 


-0.036 


-0.002 


-0.036 


-0.034 


-0.001 


0.3 


-0.065 


-0.053 


-0.018 


-0.046 


-0.040 


-0.007 


-0.040 


-0.037 


-0.003 


0.5 


-0.038 


-0.040 


0.004 


-0.044 


-0.039 


-0.006 


-0.039 


-0.037 


-0.005 


0.7 


0.089 


0.058 


0.124 


-0.023 


-0.022 


0.012 


-0.022 


-0.027 


0.008 


0.9 


0.291 


0.220 


0.331 


0.021 


0.013 


0.049 


0.013 


0.002 


0.038 



correction in the case BAREl contains both leading and next-to-leading contributions. In 
the rest cases the higher order corrections are take in the leading approximation. For a 
comparison we give in Table 2 also the corresponding numbers derived by the help of Monte 
Carlo generator BHLUMI [2]. Parameters are the same as for Table 1. The results of the non- 
exponentiated version of BHLUMI for BAREl differs from the exponentiated ones by three 
digits after decimal point, which are given in parenthesis. The numbers in square brackets are 
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Table 2 

BHLUMI results for the small-angle Bhabha cross-section (in nb) 

0{a^) beyond the first order 



Xq 


BAREl 


CALOl 


CAL02 


BAREl 


CALOl 




CAL02 




0.1 


166.046 


166.329 


131.032 


166.892 (988) [0.008] 


167.203 


131.835 


95.458 


98.834 


0.3 


164.740 


166.049 


130.739 


165.374 (471) [0.014] 


166.795 


131.450 


95.233 


98.539 


0.5 


162.241 


165.287 


130.176 


162.530 (594) [0.018] 


165.830 


130.727 


94.841 


98.020 


0.7 


155.431 


161.794 


127.528 


155.688 (620) [0.018] 


162.237 


127.969 


93.520 


96.054 


0.9 


134.390 


149.925 


117.541 


137.342 (191) [0.018] 


151.270 


118.792 


87.359 


88.554 



absolute difference (in nb) between our second order photonic correction and the one of the 
non-exponentiated BHLUMI version. BHLUMI numbers beyond the first order for CAL02 
case correspond to WW, NN and WN angular acceptances respectively; and the rest is for 
WW case. Beyond the first order all BHLUMI numbers, except the ones in parenthesis for 
BAREl, correspond to the version based on the Yennie-Prautchi-Suura exponentiation. 

On the level of the first order correction BHLUMI numbers exceed our ones approximately 
on 0.3 per mille for all variants of event selection. We think that this is due to the difference 
in our approaches: BHLUMI computes the first order correction [14] according to complete 
0{a) formulae, while we take into account only i-channel graphs as discussed above. 

To be consequent we have to compare our results due to second order photonic contribution 
with BHLUMI ones which belongs to the non-exponentiated version only. These are the 
numbers into the parenthesis for BAREl (three figure after point in the cross-section). To 
compare it needs to remove the difference due to the first order contribution. After this we 
find that our second order photonic correction which includes leading and next-to-leading 
contributions exceeds a little bit the BHLUMI one and conclude about very expressive agree- 
ment in the case of BAREl WW. 

As concerns calorimeter event selection we have not explicit calculation of the second order 
next-to-leading contribution. That is why now we can speak about the first order correction 
only. As one can see from the Tables the agreement of our numbers for WW variant of CALOl 
and CAL02 with BHLUMI ones is on the same level as for BAREl WW one. 
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